Test of multiscaling in DLA model using an off-lattice killing-free algorithm 
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We test the multiscaling issue of DLA clusters using a modified algorithm. This algorithm elim- 
inates killing the particles at the death circle. Instead, we return them to the birth circle at a 
random relative angle taken from the evaluated distribution. In addition, we use a two-level hier- 
archical memory model that allows using large steps in conjunction with an off-lattice realization 
of the model. Our algorithm still seems to stay in the framework of the original DLA model. We 
present an accurate estimate of the fractal dimensions based on the data for a hundred clusters with 
50 million particles each. We find that multiscaling cannot be ruled out. We also find that the 
fractal dimension is a weak self- averaging quantity. In addition, the fractal dimension, if calculated 
using the harmonic measure, is a nonmonotonic function of the cluster radius. We argue that the 
controversies in the data interpretation can be due to the weak self- averaging and the influence of 
intrinsic noise. 



I. INTRODUCTION 

The DLA model [l| plays the same role in the physics of 
structure growth in two-dimensions as the Ising model 
plays in the theory of phase transitions. It catches the 
main features of the random growth and is quite simple in 
definition. But this model however is still not solved ana- 
lytically more than twenty years after its introduction . 
Direct simulations of the original model and calculations 
using the conformal-mapping technique 3] are the two 
main methods for investigating DLA structures. 

It is commonly believed that DLA clusters are random 
fractals Q, Q , and the accepted estimate for the fractal 
dimension is D = 1.71 ± 0.01. The analytic result |4| 
predicts D = 17/10 in agreement with the numerical es- 
timates. The surface of a DLA cluster demonstrates the 
multifractal properties obtained in simulations 0, Q and 
supported analytically 0. 

Several groups claim DLA clusters have multiscaling 
properties |a M, EH EH [CI : that the penetration depth 
£ is scaled differently from the deposition radius i?dep 
and that a whole set of scaling exponents exists within 
the framework of multiscaling. Recently, these claims 
were doubted in papers by Somfai, Ball, Bowler, and 
Sander [H El- 

The off-lattice killing-free algorithm, our implementa- 
tion of the DLA algorithm, allows generating a large 
number of huge clusters and calculating the fractal di- 
mensions of the quantities mentioned above. Our nu- 
merical results do not support the arguments presented 
in dill but favor the early results [1 II El El EI 

The multiscaling that was "suspected" in those papers 
was attributed in 0, Ei| to the strong lattice-size effects 
they advocated, with the correction-to-scaling exponent 
1/3. We do not find any evidence for that in our data. In- 
stead, we prefer to attribute the spreading of the fractal 
dimension values as estimated from the different quanti- 
ties to the weak-self-averaging of the fractal dimension. 
Its relative fluctuation decays with the number of parti- 
cles approximately as T D = {(D 2 ) -(D) 2 ) /(D) 2 ) oc l/AT? 
with 7 w 0.35 ± 0.04 for the large cluster sizes, i.e, with 



the exponent about three times smaller than would be 
expected for the usual averaging of the quantity. Thus 
one can expect large fluctuations of the fractal dimension 
as estimated from the different quantities, with different 
methods, and from different cluster sizes. 

The alternative and more naive view is to say that 
the last exponent is also evidence for the multiscaling 
properties of DLA cluster. 

Our algorithm is off-lattice with memory organization 
similar to one used in the Ball and Brady algorithm EH • 
The main differences are: (i) we use only two layers in 
the memory hierarchy, which seems optimal for reducing 
memory usage, keeping the overall efficiency of simula- 
tions high; (ii) we use bit-mapping for the second layer 
of memory to reduce the total memory used by the sim- 
ulation program; (iii) we use large walk steps. The new 
feature is the recursive algorithm for the free zone track- 
ing. 

The essential feature of our algorithm is that we mod- 
ified the rule for the particles that go far from the clus- 
ter. We never kill any particle at the outer circle but 
return them to the birth circle pH, Il7| with an evalu- 
ated probability. This procedure eliminates the effect of 
the potential distortion (and thus of the cluster harmonic 
measure) , keeping our algorithm in the same universality 
class. 

In this paper we present all details of our algorithm. 
We believe that some details may be important both for 
understanding the results and for comparing the results 
properly. We briefly discuss the essentials of our imple- 
mentation of the off-lattice killing-free algorithm in Sec- 
tion [n] accompanied with some technical details on the 
derivation of the return-on-birth-circle probability given 
in Appendix lAl and details on the memory model in Ap- 
pendix^ Various methods for estimating the fractal di- 
mension are described in Section IIIII and compared with 
those from other off-lattice estimations. We discuss fluc- 
tuations of the fractal dimension in Section IIVI and the 
multiscaling issue in Section [V] The discussion in Sec- 
tion concludes our paper. 
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II. KILLING-FREE OFF-LATTICE 
ALGORITHM 

A cluster grows according to the following rules: (1) 
We start with the seed particle at the origin. (2) A new 
particle is born at a random point on the circle of radius 
i?b- (3) A particle moves in a random direction; the step 
length is chosen as big as possible to accelerate simula- 
tions. (4) If a particle walks out of the circle of radius 
Rd > i?b, it is returned to the birth radius i?b at the 
angle (p' taken from distribution l)A6(l and relative to the 
last particle coordinate. (5) If particle touches a cluster, 
it sticks. (6) The cluster memory is updated. Steps from 
(2) through (6) are repeated N times. 

The difference from the traditional DLA algorithm is 
in Step 4. We remind the reader that in the original 
algorithm, a particle is killed when it crosses Rd and the 
new one starts a walk from a random position on the 
circle i?b, i.e, at the position (Rh, ^random)- Clearly, the 
rule being perfect in the limit Rd — ► oo will influence the 
growth stability when Rd is finite. 

Instead, if particle crosses Rd to a position with r > 
Rd, we use the probability [lljiS as determined by ex- 
pression (|A6() in Appendix A with x — r/R^ to obtain 
new particle position. The particle then walks from the 
position (Rb,ip' + ip), assuming (r, <p) is the old particle 
position. Details on obtaining expression i|A6(l and on 
generating random numbers with given probability are 
presented in [lr3 | and in Appendix A. 

During the cluster growth, the empty space between 
cluster branches also grows. A special organization of the 
memory is implemented to avoid long walks in the empty 
space. We modify the hierarchical memory model |l5j . 
using only a two-layer hierarchy (see Appendix [5] for de- 
tails) and a bit-mapping technique for the second layer 
to reduce total memory. 



III. FRACTAL DIMENSION ESTIMATIONS 

In this section we present the results on the estimation 
for the DLA cluster fractal dimension analyzing various 
cluster lengths: deposition radius i?dcp, mean square ra- 
dius i?2, gyration radius R gyr , effective radius R c g and 
ensemble penetration depth £. The dependence of the 
length R on the number of particles N gives estimation 
of the fractal dimension through the relation R cx A 1 / 15 . 
There are two ways (see Table QJ to extract this depen- 
dence. The first is to average over the ensemble of clus- 
ters, for example, (r(N)) = Y2i=i r i{N)/K, where the 
sum is over K clusters and (N) is the position of the 
iV-th particle in the i-th aggregate. The second is to aver- 
age over the harmonic measure, which is the probability 
of sticking at the point r, for example, i?dep = / T dq. 

We use data from K = 100 clusters, each built up by 
the algorithm described in Section[n]with 5 TO 7 particles. 
A typical cluster [24| is shown in Figure ^ Estimations 
of the fractal dimension from the various cluster lengths 




FIG. 1: (Color online) Typical DLA cluster with 5 • 10 7 par- 
ticles grown using killing- free off-lattice algorithm. Color de- 
notes particle age. 
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1.70942(46) 


(Jrdq) 


1.70922(97) 


i?2 




1.71003(45) 


^Jr*dq) 


1.7087(11) 


^gyr 




1.71008(96) 






Reff 






(exp(J\nrdq)) 


1.70944(87) 




\J R2 2 — -Rdcp 2 


1.74(3) 




1.69(7) 



TABLE I: Estimates of the fractal dimension D extracted 
with the fit N oc R D to the dependence of the various lengths 
R ( deposition radius -Rdop, mean square radius R2, gyration 
radius -R gyr , effective radius 7? e ff and ensemble penetration 
depth £) on the number of particles N . The third (fifth) 
column is the fit to the data calculated with the definition 
given in the second (fourth) column. 



using both ways of averaging are shown in Tabled The 
errors given in parentheses as corrections to the last digit 
include both statistical errors and fitting errors. 

In averaging over the harmonic measure (the last col- 
umn in Table d with the preceding definition column), 
we first estimate D from a single cluster [25| and than 
average over the K = 100 samples. 

The estimations from all lengths agree well with each 
other within the error and with the most accepted value 
of the fractal dimension D = 1.711 (see, e.g. 13J). The 
error in measuring the penetration depth is much higher 
because of its complex structure. 

It was proposed in 0, Q that the various lengths 
depend on the number of particles with the correction 
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Definition 1 


Definition 2 
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Rdcp 


1.394(2) 


0.59 ± 0.28 


1.398(2) 


0.006 ± 0.030 


R 2 


1.414(2) 


0.22 ±0.27 








Rsyr 


0.964(1) 


-0.22 ± 0.01 








€ 


0.239(1) 


-13 ±1 









TABLE II: Coefficients of correction to scaling fits (expr. 
with fixed D — 1.711 and v = 0.33. The definitions of various 
length R are given in Table 



term 



R(N) = RN 1/D {1 + RN 



(1) 

and that exponent v is the same (y = 0.33) for all quan- 
tities. 

The fit of the data to expression Q with the fixed 
values of D — 1.711 and v — 0.33 is presented in TablclTl] 
Table [H] here should be compared with Table 1 in [1J| . 
The difference in the values of R is about factor of two 
and probably because of the different units of the particle 
size used in simulations. We fix to unity the particle 
radius and not the particle diameter. Result of the fit is 
extremely sensitive to the value of D used. For example, 
if we fix fractal dimension D to the value 1.710 which 
one may suggest from our Table |I| the values of R and R 
for the fit to R gyr changed from those in the third line 
of our Table|ll|to 0.958(1) and -0.06(1). Thus, values of 
R differ by six standard deviations, and values of R - by 
sixteen standard deviations. 

If we fit our data for the different cluster realizations 
to the expression (JIJ without fixing D and u, then we 
find a large fluctuation of R around the zero value. 

If we suppose that authors of Ref. fixed diameter to 
unity, we may conclude that our data for R from Table ITT| 
coincide with the corresponding data in Table of Ref. [TJ( 
and not the data for R. We can therefore conclude that 
the results of fitting to expression |QJ are inconclusive and 
that the values of the coefficient R presented in Table |n] 
are just random. 

These large fluctuations can be understood in the 
framework of the weak self-averaging of the fractal di- 
mension, which we describe in the next section. 



IV. WEAK SELF-AVERAGING OF D 

We check how the fluctuations of the measured frac- 
tal dimension D depend on the system size N. By 
analogy with thermodynamics, the relative fluctuation 
T D = '{D 2 ) - {D) 2 )/{D) 2 of the quantity D should de- 
crease as the inverse system size. For full self-averaging 
of the fractal dimension, Tu oc 1 /N oc 1/ R D can be ex- 
pected. In the case of a slower decay Td <x 1/N 1 with 
7 < 1, one can say that self- averaging of D still occurs 
and that it is weak self-averaging. 



We extract the fractal dimension from the analysis of 
the clusters in two ways. First, we calculate the number 
N of particles inside the circle of radius R for the given 
cluster. The slope of this curve on a log-log plot gives Di. 
The fractal dimension Di as a function of TV is denoted as 
Di(N), where i is the number of cluster, i = 1,2, K, 
and K = 100. Then we average Di(N) over the ensem- 
ble of K clusters, D(N) = ± J^iLi A (AO- Th e fractal 
dimension D(N) is plotted in Figure [21 with a bold line 
as a function of the system size N. To better under- 
stand the behavior of D, we also present the results of 
averaging over smaller ensembles. We divide the whole 
ensemble with 100 clusters into five independent groups. 
Averaging over each group gives five different curves for 
D(N), which exhibit strong fluctuations. Error bars are 
computed as fluctuations of Di(N) in the ensemble of 
K = 100 clusters. 

For sufficiently large N > 10 5 , the values of D{N) 
vary mainly in the range 1.695-1.715, which is about the 
usually accepted value of the fractal dimension. The drop 
off of the curve D(N) at N > 2-10 7 is due to the influence 
of the cluster boundary: the most active zone of cluster 
growth, which is underdeveloped in comparison with the 
rest of the cluster, is now inside R. 

Relative fluctuations Td of that quantity are shown in 
Figure El The bold line represents 100-cluster ensemble 
averaging compared with five lines computed using five 
20-cluster groups. Fluctuations decrease with the expo- 
nent 7 = 0.33±0.02. This is three times slower decreasing 
than expected in the case of full self-averaging. 

Next, we analyze self-averaging of the fractal dimen- 
sion as extracted from the dependence of the deposition 
radius Rde P (N) calculated by averaging over the har- 
monic measure. This quantity averaged over the ensem- 
ble of 100 clusters is plotted in Figure ||as a function 
of the system size -Dh a rm(-/V). It exhibits some "oscilla- 
tion" around the value 1.708, which is very close to the 
accepted DLA fractal dimension. Averages over smaller 
ensembles also demonstrate this feature. It is not clear 
how this quantity will change as the system size increases 
further. 

Relative fluctuations ^r; harm of Z?harm are shown in 
Figure |5] There are two regimes of the Jc hllm decay. 
First, for the cluster sizes N < 10 6 , it decays with 
7 = 0.71 ± 0.02, much faster than for the To as esti- 
mated with the conventional counting method described 
above with 7 w 0.33. The next regime, for the larger 
system sizes N > 10 6 , shows slower decay with the ex- 
ponent 7 = 0.38 ± 0.01, close to those estimated by the 
traditional counting method 26] . 

In practice, the exponent value 7 ~ 0.33 means that to 
obtain more accurate estimate of the fractal dimension, 
one have to increase the number of particles in the cluster 
by three orders of magnitude. 
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FIG. 2: (Color online) Fractal dimension D(N) as a function 
of the number TV of particles inside the radius R. Solid curve 
with error bars represents D{x) averaged over 100 clusters. 
Another curves are the averages over five different 20-clusters 
ensembles. 
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FIG. 3: (Color online) Decay of the relative fluctuations of 
the fractal dimension shown in Figurc|5| The symbols are the 
same as in Figure |5] The solid line in the inset is the linear 
fit to the open circles. 



V. MULTISCALING 

It was l ong di scussed that DLA clusters are not simple 
fractals [8ll9l ll0llrTllT^ | and, for example, that the fractal 
dimension depends continuously on the normalized dis- 
tance from the cluster origin It was suggested in [llj 
that the density of particles at a distance r from the ori- 
gin obeys the equation g(r, R gyr ) = c (x)R^y^ d ', where 
x = r I i? gyr , with a nontrivial (non constant) multiscaling 
exponent D (x) (lines with open symbols in Figure 0. 

Quite recently, Somfai et al. [l3L FLU ] claimed that the 
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FIG. 4: (Color online) Fractal dimension .Dharm(TV) as a func- 
tion of the number of particles in the cluster used to estimate 
Rdep with averaging over the harmonic measure. The symbols 
are the same as in Figure [5] 
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FIG. 5: (Color online) Decay of the relative fluctuations of the 
fractal dimension .Dharm(TV) shown in Figure 01 The symbols 
are the same as in Figure [5] Inset: the linear fit to the data 
with TV up to 5 ■ 10 5 and the linear fit to the data with larger 
TV. 

multiscaling picture is wrong and is misled by finite-size 
transients. They argued for a strong dependence of the 
radius estimators (Rdcp,R2, etc., see Table ITT)l on the 
system size Expr. l|T)l. where the leading subdominant 
exponent is estimated as v « 0.33. 

It is well known 0, 0, |^| that the fractal dimen- 
sion D(x) can be found using the probability P(r, TV) 
for the N-th particle to be deposited within a shell of 
width dr at a distance r from the seed. The simplest 
and most obvious form of the probability is P(r, TV) = 

WT)f( ~^Wr~^ )) wnere fiv) i s Gaussian. Practically, 
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the Gaussian distribution can be obtained by averaging 
over a large number of clusters. For the single-cluster re- 
alization this function has some particular form reflecting 
the details of the cluster growth as shown in Figure [3 in 
which each local maximum is associated with an actively 
growing branch. 

In [14| Somfai et al. computed D(x) from the Gaussian 
probability P(r, TV), and D(x) with corrections to scaling 
coincides well with the numerical results of Amitrano et 
al. 0| • Somfai et all state that D(x) tends to a constant 
value as TV — > oo. In other words, Somfai et al. advocate 
that multiscaling is transient and is an artifact of the 
finite size of the DLA clusters. 

In contrast, our numerical results demonstrate that 
(1) there is no evidence for the finite-size corrections with 
the exponent v = 0.33, (2) D(x) seems not tend to a con- 
stant, and (3) it is not correct to use the Gaussian prob- 
ability P(r, TV) to compute D(x) for the DLA model. 

Gaussian distribution does not reflect details of the 
DLA cluster because it is the outcome of the the aver- 
aging over a large number of clusters. After such aver- 
aging all details of the random nature of the growth of 
a particular cluster are washed out, and the Gaussian 
distribution is just the result of the central limit theo- 
rem. There are some indications by Hastings Q, who 
computed DLA fractal dimension from field theory, that 
the Gaussian distribution by itself is insufficient for de- 
scribing DLA clusters, and some noise must be added 
to obtain a model corresponding to the DLA model. To 
some extent, the local maxima in Figure [7] are due to 
that noise, in contrast to the distribution averaged over 
the cluster ensemble, which is smooth. 

Accordingly, D(x) calculated from the averaged distri- 
bution would be constant in the limit of large TV. We can 
therefore say that multiscaling comes from the fluctua- 
tions (natural noise) of the DLA cluster growth process. 

The fractal dimension D(x) for different cluster sizes 
are shown in Figure|S]together with the results from [Tl) . 
There is a notable maximum in D(x) at around x = 1.4, 
and its size does not change significantly with the number 
of particles in a cluster. We note that the position of the 
maximum can be found from the ratio of (i?dcp — £) to 
i? gyr . In our simulations, it equals 1.19, which coincides 
well with the data in Figure 

The maximum seems to occur because growth is mostly 
completed in the region r < i?d cp — £ , and the fractal di- 
mension D(x) decrease to the left of x = (Rd cp — £)/i? gyr 
in the actively growing region D(x) because of the lack 
of particles there. For small x, i.e., for r << i? gyr , no 
new particles are being added, and D(x) is the same for 
different cluster sizes. 

We note that the prediction of Somfai et al. for D(x) 
with a system size TV = 10 7 is quite smaller than the one 
computed by us and presented in Figure We also es- 
timate limit of the three curves plotted in Fig. H3 taking 
limit of iV — ► oo at the fixed value of x, and plot result 
with the solid circles in Fig. Thus, our data does not 
demonstrate tendency of D(x) to a constant value, and 




FIG. 6: (Color online) Multiscaling fractal dimension D(x) for 
different cluster sizes: solid triangles, TV = 10 6 ; solid squares, 
TV = 10 7 ; stars, N — 5 ■ 10 7 ; solid circles, limit of our data for 
TV — » oo; open symbols are from [TH| : open squares, TV = 10 4 
square lattice; open circles, TV = 10 5 off-lattice. 



rather support our observation that the dynamic growth 
of the cluster is dominated by the active zone, and max- 
imum in D(x) reflects this nature of DLA. At the same 
time we have to note that accuracy of the data presented 
here is not enough for the final decision, and future inves- 
tigation with the higher accuracy have to be still done. 



VI. SUMMARY 

We have implemented modifications to the DLA algo- 
rithm that help us to grow large DLA clusters and test 
some recent claims about its properties. In our exper- 
iments, aggregates do not exhibit corrections to scaling 
laws. Nevertheless, the results show multiscaling prop- 
erties. This means that there should be another way for 
such a clusters to appear. We tried to analyze the pro- 
cesses responsible for multiscaling. We will address this 
question in future research. 
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FIG. 7: Probability P(r) for a particle to stick at the distance 
r, computed for a single DLA cluster. 



APPENDIX A: PROBABILITY TO BE ALIVE ON 
THE BIRTH CIRCLE 

The main result of this Appendix is expression HA6II. 
The same expression was obtained earlier [27| in Refs. [la, 
Il7j . We found our result is still worth publishing, since 
we derive it from a different point of view. 

We consider a particle at some position (r; <f> = 0) out- 
side circle of radius Rb, T > Rb- The particle moves 
randomly, and its size and step is much smaller than r 
and i?b • The question is what is the probability for a par- 
ticle to intersect the circle Rb at the angle ip' . Clearly, 
for the particles walking from infinity, i.e., r >> Rb, that 
distribution is uniform in [— ir;ir]: 



P{f)d<p = -^-dtp. 

ITT 



The conformal map 



w = R] 



zr + R h 
' r + zRb 



(Al) 



(A2) 



maps infinity to (r, 0) and the unit circle to the circle 
of radius Rb- Transformation (|A2|) changes probabil- 
ity IjAlfl to the modulus of derivative of the conformal 
map, 



dz 
dw 



Rb{r 



w 



(A3) 



Substituting w = Rb exp(i<//) in I) A3|l . we obtain the re- 
sulting probability 



const 



x 2 — 2x cos ip' + 1 



(A4) 



as a function of the ratio x = r/Rb > 1. This is a 
probability for the particle beginning its walk at the point 



FIG. 8: Sketch of probability P(ip') of the walk intersection of 
the circle Tib with the initial position of the particle at (r, 0). 



(r, 0) outside Rb to intersect circle Rb at the angle ip' (see 
Figure IHJ. 

The constant in (|A4(I is associated with the probability 
for the particle to move to infinity. It can be identified 
using the analogy between the DLA model and Laplacian 
growth as was pointed in |l8[ for the dielectric breakdown 
model, which is a generalization of the DLA model. 

The probability that the particle sticks somewhere to 
the cluster is proportional to the electrical field at that 
point. We consider three circles of radii i?2, and 
i?3, i?i < i?2 < ^?3- The external circles R\ and R3 are 
set under the potential $ = 0, and particles stick one 
of them. The circle R2 (place of a birth for particles) 
is set with $ = 1. The solution of the Laplace equa- 
tion V 2 $ = with the above boundary conditions gives 
the distribution of the potential and thus of the electri- 
cal field E. The probability for the particle to stick on 
the circle R\ would be P\ = const * /i r i_ Rl |E|d 2 r and 
similarly for P3. The ratio of these probabilities is: 



P3 
Pi 



logf 



(A5) 



The probability P3 vanishes as i?3 ^ 00: P3 — > 0. 
This means that all particles starting on the birth ra- 
dius should irreversibly collide with a cluster and the 
constant in (|A4|) is easily found by normalizing the prob- 
ability, const = (x 2 — 1)/2tt. The final expression for the 
probability 



1 



x 2 - 1 



27T x 2 — 2x cos ip' + 1 



(A6) 



provides correct values for both the limits x — > 00 and 
x — v 1. The same expression for probability was de- 
scribed in , but it was obtained other way. Transfor- 
mation 



f(u) = 2 arctan 



1 



tan(wr/2) 
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maps uniformly distributed in interval [—1,1] random 
number u to random number with distribution (|A6() . 



APPENDIX B: MEMORY ORGANIZATION 

Step 3 of the algorithm described in Section [H] es- 
sentially contains two routines. We must choose, first, 
the direction of the random walk and, second, the walk 
length. The direction of the walk is chosen uniformly 
in [0;27r]. Because the motion is uniform in direction, 
we can increase the step, but only if particle is far away 
from the cluster (in our simulations we choose that dis- 
tance such that particle is more than five units away from 
the cluster, otherwise, its step is one unit of length). The 
length for the big step is chosen with the condition that 
the walk should not intersect any particles of the cluster. 
Therefore, the distance to the cluster d pc is evaluated, 
and the step length is taken as that distance. 

The reason for implementing variable step length is 
as follows. Most of the time the particle is moving far- 
away from the cluster, and choosing a step length of the 
order of the particle distance from the cluster accelerates 
simulations. 

To realize all the proposed improvements efficiently we 
must organize memory in a special manner. When the 
particle moves, we must know whether it collides with a 
cluster, and to check this, we must iterate over all par- 
ticles. Such an approach is rather unreasonable and we 
must therefore restrict the number of particles to test. 
This is easily done by dividing the space into square cells, 
each about twenty units of length. Each cell saves infor- 
mation about particles that stuck in the region covered 
by it. We therefore need only check cells that are in the 
region of one step. 

This model also improves the process of seeking the 
size of free space for the big step. To find the distance 
from current position to the cluster precisely is a difficult 
task, but it suffices in most cases to know it with an accu- 
racy estimated from below, the size of one cell. Figure 
shows how this is done. Cells are plotted in the picture as 
squares with bold lines. They are marked with numbers 
showing their distance from the particle location. For 
simplicity, cells with the same number are thought to be 
at the same distance from the particle. Occupied cells are 
shaded. In this example the particle is allowed to jump 
with the step length L — 2R, where L is the cell size, and 
R the particle radius. The distance L is the radius of the 
inscribed circle with a center somewhere in cell (in the 
worst case, the center lies on the border of cell 0). This 
length should be reduced because cells save only particle 
centers and there could be a projection of a particle into 
another cell with a size of R. To simplify the algorithm 
we seek the step length only using the free/occupied cells 
picture. The first step is to check whether there are other 
particles in the cell we are now in, then we should check 
cells marked with 1, then marked with 2 (not shown in 
picture), and so on until we find an occupied cell. 




FIG. 9: Determination of the step length. 



During DLA growth, the intervals between branches 
increase notably, and the time to traverse all cells while 
seeking an occupied one also increases. To reduce the 
influence of such a process, Ball and Brady 0] developed 
a hierarchical memory model, where one cell is divided 
into smaller subcells and so on. This approach seems 
rather memory consumptive: in growing a large cluster, 
it would become a bottleneck of an algorithm. 

The desired effect can be achieved in another way. Be- 
cause the cell size is much more bigger than the par- 
ticle size, the distribution of free cells changes slowly, 
and process of seeking the maximum free space could 
be started not from the particle position but from the 
free line achieved in the previous search from this origin. 
To implement this we must save the value of free space 
around each cell. If this information is unknown, i.e., it is 
the first time to seek the maximum step length from the 
current cell, we should traverse all cells from the begin- 
ning; otherwise we start from the line previously saved 
and move to the center. 

The cell size is chosen as follows. It should not be very 
small: a small size results in high memory consumption 
and a rapid change of the distribution of free cells. On the 
other hand, its size restricts the precision of the length 
we find for the big step, i.e., there is a region of cell size 
near the cluster where particles can move only with a 
small steps. To avoid this last constraint, we use second- 
layer information. Like Ball and Brady, we divide the 
cells into smaller ones. To minimize memory consump- 
tion, we realize them as an 32 bit integer, where each 
bit shows whether the corresponding subcell is occupied. 
Each big cell can therefore be divided into not more than 
25 subcells. 
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As mentioned above, the second-layer information is 
only used when the particle moves near the cluster, where 
the accuracy given by the first layer is insufficient. Fig- 
ure shows large cells (bold borders) and small cells. 



Using the first layer we can find that only the space in- 
side the dashed circle is free. The second layer gives a 
more precise result: the particle can jump up to the bold 
circle. 



[1] T.A. Witten and L.M. Sander, Phys. Rev. Lett. 47 (1981) 
1400. 

[2] A. Bunde and S. Havlin, eds., Fractals and Disordered 

Systems (Springer, Berlin, 1996). 
[3] M.B. Hastings and L.S. Levitov, Physica D 116 (1998) 

244. 

[4] M.B. Hastings, Phys. Rev. E 55 (1997) 135. 

[5] P. Meakin, H.E. Stanley, A. Coniglio, T.A. Witten, Phys. 
Rev. A 32 (1985) 2364. 

[6] T.C. Halsey, P. Meakin, and I. Procaccia, Phys. Rev. 
Lett. 56 (1986) 854. 

[7] T.C. Halsey, B. Duplantier, and K. Honda, Phys. Rev. 
Lett. 78 (1997) 1719. 

[8] M. Plischke and Z. Racz, Phys. Rev. Lett. 53 (1984) 415. 

[9] B. Davidovitch, H.G.E. Hentschel, Z. Olami, I. Procac- 
cia, L.M. Sander, and E. Somfai, Phys. Rev. E 59 (1999) 
1368. 

[10] A. Coniglio and M. Zannetti, Physica A 163 (1990) 325. 
[11] C. Amitrano, A. Coniglio, P. Meakin, and M. Zannetti, 

Phys. Rev. B 44 (1991) 4974. 
[12] P. Ossadnik, Physica A 176 (1991) 454, ibid. 195 (1993) 

319. 

[13] R.C. Ball, N.E. Bowler, L.M. Sander, E. Somfai, Phys. 

Rev. E 66 (2002) 026109. 
[14] E. Somfai, R.C. Ball, N.E. Bowler, L.M. Sander, Physica 

A 325 (2003) 19. 
[15] R.C. Ball and R.M. Brady, J. Phys. A 18 (1985) L809. 
[16] E. Sander, L. M. Sander, R. M. Ziff, Computers in 

Physics 8 (1994), 420-425; L.M. Sander, Contemporary 



Physics 41 (2000), 203-218. 

[17] H. Kaufman, A. Vespignani, B.B. Mandelbrot, and L. 
Woog, Phys. Rev. E 52 (1995) 5602. 

[18] L. Niemeyer, L. Pietronero, H.J. Wiesmann, Phys. Rev. 
Lett., 52 (1984) 1033. 

[19] P. Meakin, J. Phys. A 18 (1985) L661. 

[20] P. Meakin, L.M. Sander, Phys. Rev. Lett. 54 (1985) 2053. 

[21] J. Lee, S. Schwarzer, A. Coniglio, H.E. Stanley, Phys. 
Rev. E 48 (1993) 1305. 

[22] S. Tolman and P. Meakin, Physica A 158 (1989) 801; 
Phys. Rev. A 40 (1989) 428. 

[23] T.A. Rostunov and L.N. Shchur, JETP 95 (2002) 145. 

[24] a high quality figur e is available at 

http: / /www. comphys.ru / may /dla-cluster .eps 

[25] Averaging over the harmonic measure in practice is done 
by freezing the growth and using probe particles that 
walk according to the usual growth algorithm but do not 
stick to the cluster. Then a new particle starts its mo- 
tion. We use 10 4 probe particles typically. The average 
of the positions r where they touch the cluster gives the 
quantity averaged over the harmonic measure. 

[26] We note that the close value of the exponent 7 was es- 
timated for the quasi-flat growth of a DLA cluster using 
Hastings-Levitov conformal mapping simulations |23|. 

[27] We thank R. Ziff and anonymous referee for pointing us 
this papers. Authors of the present paper were not aware 
on that result prior the paper submission. 



